Signature of a three-dimensional photonic band gap observed on sihcon inverse 

woodpile photonic crystals 
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We have studied the reflectivity of CMOS-compatible three-dimensional silicon inverse woodpile 
photonic crystals at near-infrared frequencies. Polarization-resolved reflectivity spectra were ob- 
tained from two orthogonal crystal surfaces using an objective with a high numerical aperture. The 
spectra reveal broad peaks with maximum reflectivity of 67% that are independent of the spatial 
position on the crystals. The spectrally overlapping reflectivity peaks for all directions and polar- 
izations form the signature of a broad photonic band gap with a relative bandwidth up to 16%. 
This signature is supported with stopgaps in plane wave bandstructure calculations and with the 
frequency region of the expected band gap. 



I. INTRODUCTION 

Currently, many efforts are devoted to create an in- 
tricate class of three-dimensional meta-materials known 
as photonic crystals that radically control propagation 
and emission of light [Ilfl2]. Photonic crystals are or- 
dered composite materials with a spatially varying di- 
electric constant that has a periodicity of the order of 
the wavelength of light [H] . As a result of the long-range 
periodic order, the photon dispersion relations are orga- 
nized in bands, analogous to electron bands in solids 
Frequency ranges called stopgaps emerge in which light 
is forbidden to propagate in particular directions due 
to Bragg interference [TS]. In specific three-dimensional 
crystals, a common stopgap is formed for all polarizations 
and directions, called the photonic band gap. Light can- 
not propagate inside the photonic band gap, allowing for 
ultimate control of light emission. Emission rates and di- 
rections can be manipulated [5H51[T^. which could lead to 
efficient micro-scale light sources O [H] and solar cells 
[TT] . Additional interest is aroused by the possibility of 
Anderson localization of light by point defects added to 
photonic band gap crystals P]. 

It is an outstanding challenge to experimentally 
demonstrate a photonic band gap; inside the photonic 
band gap the density of optical states equals zero. The 
density of states can be investigated with light emitters 
placed inside the crystal, see e.g. Rcf. [6-8J. These ex- 
periments are difficult to perform and to interpret, and 
are not always possible due to the limited availability of 
appropriate light sources and detection methods. How- 
ever, one can obtain an indication of the band gap by 
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observing a stopband in a directional experiment, such 
as a peak in reflectivity or a trough in transmission |181 - 
I23j . The frequency width of an experimentally observed 
stopband often corresponds to a stopgap in the disper- 
sion relation, where light is forbidden to propagate. In- 
terestingly however, a peak in reflectivity or a trough 
in transmission also occurs when an incident wave can- 
not couple to a field mode in the crystal [131 HH l25] . 
Therefore, one typically compares observed stopbands 
with theory such as calculated bandstructures to indi- 
cate the presence of a band gap. Unfortunately, most 
theory is typically valid for ideal structures, for instance 
of infinite size. More compelling evidence of a photonic 
band gap could be obtained if one demonstrates that 
stopbands have a common overlap range independent 
of the measurement direction, since light is forbidden 
to propagate in a photonic band gap. Angle-resolved 
PUI US] and angle-averaged spectra [23J have been col- 
lected, however, due to experimental considerations it 
appears to be extremely difficult to measure over At: sr 
solid angle. Stopbands are typically investigated on only 
one surface of the photonic crystal, and only few studies 
have investigated multiple surfaces, see e.g. Ref. |27] . In 
addition, one needs to rule out spurious boundary effects 
by confirming that stopbands reproduce at different lo- 
cations on the crystal, requiring position-dependent ex- 
periments. Furthermore, it can occur that field modes 
in the crystal can only couple to a specific polarization 
[28] . Polarization-resolved experiments are required to 
demonstrate that stopbands are present for all polariza- 
tions [29'. Therefore, a strong experimental signature for 
a photonic band gap is obtained if one can demonstrate 
that stopbands are position-independent and overlap for 
different directions and orthogonal polarizations. To the 
best of our knowledge, such a detailed analysis of stop- 
bands for three-dimensional photonic band gap crystals 
has not yet been reported. 
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In this paper we study silicon three-dimensional in- 
verse woodpile photonic crystals [30l |3T]. The inverse 
woodpile photonic crystal is a very interesting nanopho- 
tonic structure on account of its broad theoretical pho- 
tonic band gap with more than 25% relative gap width. 
Schilling et al. [23] were the first to investigate stop- 
bands in inverse woodpile crystals. They measured un- 
polarized reflectivity along one crystal direction using an 
objective with numerical aperture NA = 0.57. An indica- 
tion for the photonic band gap was found by a stopband 
that agreed with the calculated band structure. Other 
groups have fabricated inverse woodpile photonic crys- 
tals of different materials using several methods, and also 
performed reflectivity measurements along one direction 
with unpolarized light |34H36| . resulting in similar indi- 
cations of the band gap. Here, we present an extensive 
set of polarization-resolved reflectivity spectra of silicon 
inverse woodpile photonic crystals. We have collected 
spectra from two orthogonal crystal surfaces using an ob- 
jective with a high numerical aperture NA = 0.65. The 
reflectivity spectra were obtained on different locations 
on the photonic crystal surfaces to confirm the repro- 
ducibility, to determine the optical size of the crystals, 
and to investigate boundary effects. We demonstrate 
position-independent overlapping stopbands for orthog- 
onal polarizations and crystal directions, which is a sig- 
nature of a three-dimensional photonic band gap. This 
signature agrees with calculated stopgaps in plane wave 
bandstructure calculations and with the frequency region 
of the expected band gap. 



II. INVERSE WOODPILE PHOTONIC BAND 
GAP CRYSTALS 

Fig. ^a) illustrates the orthorhombic unit cell of an 
inverse woodpile photonic crystal (left), together with a 
crystal consisting of eight unit cells (right). The structure 
can be described by two identical sets of pores with radius 
r running in two orthogonal sides of a box, where each 
set represents a centred-rectangular lattice with sides a 
and c. If the ratio ^ equals ^/2, the crystal symmetry 
is cubic and the crystal structure is diamond-like, see 
Ref. [SDHSl]. This type of inverse woodpile photonic 
crystals with | = \/2 has a maximum band gap width 
of 25.4% for ^ = 0.245. The sets of pores are oriented 
perpendicular to each other and the centrers of one set 
of pores are aligned exactly between columns of pores of 
the other set, resulting in the structure of Fig. [ija). In 
fig. ^a) a coordinate system is introduced that is used 
in this paper. 

We have fabricated multiple inverse woodpile photonic 
crystals with different pore radii in monocrystalline sili- 
con with a CMOS compatible method, as is described in 
detail in Ref. [38l[39]. Fig. ^b) shows a scanning elec- 
tron microscope (SEM) image of one of these crystals 
from the same perspective as in Fig. [ija). This crystal 
consists of more than 10'^ unit cells and is located on top 






Figure 1: (color online) Structure of inverse woodpile pho- 
tonic crystals. (a). Schematic representation of the or- 
thorhombic unit cell of a cubic inverse woodpile photonic 
crystal (left) and a crystal consisting of eight unit cells (right). 
The structure can be described by two identical sets of pores 
running in two orthogonal sides of a box. (6), Scanning elec- 
tron microscope image of a cubic inverse woodpile crystal with 
a = 693 ± 10 nm, c = 488 ± 11 nm (f = ^) and r = 145 ± 9 
imi, surrounded by a macroporous two-dimensional crystal. 
The dashed orange line indicates a part of the boundary of 
the inverse woodpile photonic crystal, (c). Scanning electron 
microscope image viewed from the FX-direction. The rectan- 
gle represents one face of the orthorhombic unit cell, compare 
with (a). 



of bulk silicon. The crystal is surrounded by a macro- 
porous two-dimensional photonic crystal that forms the 
first set of pores used to fabricate inverse woodpile pho- 
tonic crystals. A second set of pores is oriented perpen- 
dicular to this two-dimensional crystal to form the inverse 
woodpile photonic crystal. A part of the boundary of the 
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inverse woodpile photonic crystal is marked with the or- 
ange dashed line. The inverse woodpile photonic crystal 
has lattice parameters a — 693 ± 10 nm, c = 488 ± 11 nm, 
and a pore radius of r = 145 ± 9 nm in both directions 
[37j . Typically eight different inverse woodpile photonic 
crystals are made on one two-dimensional photonic crys- 
tal that each extend over approximately 7 /xm x 5 /^m x 
5 ^m m size. Fig. [ijc) shows a SEM image of a crystal 
viewed parallel with the second set of pores or parallel 
with the rX-direction. In Ref. [38j it has been deter- 
mined that the second set of pores are precisely centred 
to within Ay = 17 ± 12 nm between columns of pores of 
the first etch direction, which is extremely close to the 
ideal structure of Fig. ^a) . One unit cell is marked in the 
right bottom corner of Fig. [ijc). Here, we describe an 
extensive set of reflectivity measurements collected from 
this inverse woodpile photonic crystal centred along the 
TX- and FZ-direction, which is representative for more 
than five other crystals that we have studied. 

Fig. [2] shows the bandstructure for light in a silicon 
inverse woodpile photonic crystal with ^ = 0.190 and 
tsi = 12.1 for the irreducible Brillouin zone of a sim- 
ple orthorhombic lattice, calculated with the method of 
Ref. This bandstructure appears to be representa- 

tive for the crystal of Fig [l|a) . Between reduced fre- 
quency f = 0.395 (5727 cm"!) and f = 0.460 (6668 
cm~^) a broad photonic band gap appears with 15.1% 
relative width (blue bar). The bandstructure in both the 
TX- and FZ-direction is identical and shows two stop- 
gaps, marked by the grey rectangles in the FZ-direction. 
The lowest frequency stopgap between ^ = 0.310 (4505 
cm~^) and ^ = 0.318 (4611 cm~^) is narrow and it is 
closed when going in the Z[/-direction. The broad sec- 
ond stopgap between reduced frequency ^ = 0.395 (5727 
cm-i) and f = 0.488 (7062 cm"!) is part of the pho- 
tonic band gap. The low frequency edge of the photonic 
band gap is bounded along the TX- and FZ-direction, 
which we observe in our experiment. The high frequency 
edge of the band gap occurs at the S and T points, which 
are not accessible in our experiment. Deviations in the 
crystal geometry affect the dispersion relation and there- 
fore the band gap [311 [32] . The band gap of this kind of 
structure is robust to most types of deviations, whose tol- 
erances are well within reach of our fabrication methods 



III. EXPERIMENTAL SETUP 

A supercontinuum white light source (Fianium SC- 
450-2) with a frequency range of 4000 to 22000 cm^^ 
is used to illuminate the photonic crystals. Light is po- 
larized and focused on the crystal using a gold-coated re- 
flecting objective to avoid dispersion (Ealing X74). The 
numerical aperture NA=0.65 of the objective results in a 
spectrum angle-averaged over approximately 0.447r±10% 
sr solid angle in air. The diameter of the focused beam 
is estimated to be 2wq « 1 /im from experiments on 
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Figure 2: (color online) Bandstructure of an inverse woodpile 
photonic crystal. The dispersion relation for the eight lowest 
frequency Bloch modes are calculated for a structure with a = 
690 nm, c = 488 nm and pore radius r = 131 nm (^ = 0.190). 
The left axis represents the absolute frequency in cm~^, the 
right axis shows the reduced frequency. The blue bar marks a 
photonic band gap with a relative gap width of ^ — 15.1%. 
The grey bars mark stopgaps in the FZ-direction. 



micropillars |33j . Reflectivity is measured for a broad 
range of wavevectors that are centred on the TX- and 
FZ-directions, which is expected to result in similar stop- 
bands because of the symmetry of the crystal. Reflected 
light is collected by the same objective, and the polariza- 
tion is analysed. The spectrum is resolved using Fourier 
transform infrared spectroscopy (BioRad FTS-6000) in 
combination with an external InAs photodiode. The re- 
flectivity of crystals is collected between 4000 and 10000 
cm~^ with a resolution of 32 cm~^. For calibration, 
spectra are normalized to the reflectance spectra of a 
gold mirror, which were collected before and after mea- 
suring on a crystal. The individual gold spectra show 
only minor differences due to the excellent stability of 
the setup, except between 9200 — 9600 cm^^, close to 
the master frequency of the white light source, which is 
thus excluded. The reflectivity is measured for two or- 
thogonal polarizations, where the orientation of the set 
of pores perpendicular to the measurement direction is 
used as reference; light is _L-polarized when the electric 
field is perpendicular to the direction of this set, light 
is 1 1 -polarized when it is parallel with the direction of 
this set. For example, when the reflectivity is measured 
centred on the FAT-direction (see Fig. [ija)), light is ||- 
polarized when the electric field is parallel with the FZ- 
direction. The samples are placed on a three-dimensional 
x-y-z-translation stage (precision ±50 nm) to study the 
position-dependence of the reflectivity. 



IV. RESULTS 

Position-dependent reflectivity experiments have been 
performed by scanning the photonic crystal through the 
focal spot. These scans were performed on two orthog- 
onal surfaces of the crystal allowing for measurements 
along the FA'- and FZ-direction, see Fig. [ijfo). These 
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scans provide insight in the reproducibiUty of the re- 
flectivity, the optical size of the crystal, and boundary 
effects. We will first describe one specific scan for the 
reflectivity obtained centred on the FX-direction for _L- 
polarized light. Next we will compare reflectivity spectra 
obtained from different scans to demonstrate position- 
independent overlapping stopbands for orthogonal polar- 
izations and crystal directions that are the signature for 
the photonic band gap. 



A. Position-dependent reflectivity along the 
rX-direction for _L-polarized light. 

Fig. [3] shows the position-dependent reflectivity along 
the FX-direction for _L-polarized light, where the focal 
spot was moved in the y-direction. Fig. [sj^a) shows an 
optical microscope image of the inverse woodpile pho- 
tonic crystal (dark square in the middle) taken during 
the reflectivity scan, viewed along the same perspective 
as Fig. [TJc). The bright spot centred on the inverse 
woodpile photonic crystal is the focal spot of the super- 
continuum white light source. The scan was started with 
the focal spot on the surrounding two-dimensional crys- 
tal to the left of the inverse woodpile photonic crystal, 
continued over the inverse woodpile photonic crystal and 
finished on the two-dimensional crystal to right of the 
inverse woodpile photonic crystal. Fig. ^b) shows a 
contour-plot of the measured reflectivity as a function 
of position. Fig. [sf^c) shows three cross-sections of the 
contour-plot of Fig. ^b). 

In Fig. ^b) a broad reflectivity peak is observed for 
the inverse woodpile photonic crystal {6 < y < 14 /xm) 
between 5600 and 6900 cm~^ with a high reflectivity of 
60%. The spectrum taken at y = 11 /im is shown in 
Fig. [Sf^c). This distinct peak is observed in the same 
frequency range over the entire three-dimensional crystal 
surface. The stopband corresponds to the stopgap in the 
rX-direction, which is marked by the blue bar in Fig. 
[sjjc), and which is part of the photonic band gap. 

The maximum reflectivity is about 60% for our inverse 
woodpile photonic crystals, which is likely limited by sur- 
face roughness that scatters light in directions that are 
not collected by the objective. This assertion is sup- 
ported by the observation of lower reflectivity on crystal 
surfaces that were less well polished. The finite thick- 
ness of the crystals could also lead to some reduction 
of the reflectivity. From the full width at half maxi- 
mum of the peak, we derive a Bragg length equal to 
approximately 7d, with d the spacing between lattice 
planes [H]. Therefore, we conclude that our crystals 
have a thickness of approximately 1.5 times the Bragg 
length. This corresponds to a maximum reflectivity of 
1 — e~^'^ — 78%. Since the observed maximum reflec- 
tivity is consistent with surface roughness and finite size 
effects, it is reasonable that the crystal has a high re- 
flectivity through all of the OAiir ± 10% sr probed solid 
angle. As is seen in the bandstructures in Fig. [2] the 



stopgap hardly shifts in frequency for different directions 
in the Brillouin zone, hence it is even more likely that 
the obtained angle-averaged reflectivity is representative 
for the full probed range. 

From the scan in Fig. ^b), we observe a clear dis- 
tinction between the reflectivity of the three-dimensional 
inverse woodpile photonic crystal and the surrounding 
two-dimensional crystal. A reflectivity maximum of 70% 
is observed on the two-dimensional crystal at y < 6 /im 
and at 2/ > 14 fim. The spectrum taken at y — 0.42 /im 
is shown in Fig. ^c). This peak belongs to the TK stop- 
gap between 4700 and 6300 cm~^ for TE polarized light 
[32]. At the interfaces between the two-dimensional and 
inverse woodpile photonic crystal at y = 6 and at y = 14 
/im, a peculiar double peak appears as shown in Fig. 
[sj^c), which is only present at the two-dimensional crys- 
tal to three-dimensional crystal interface. The spectrum 
at y = 6 /im is shown in Fig. [sf^c). This double peak 
consist of one component between approximately 4700 
and 5500 cm~^, and a second component between ap- 
proximately 5500 and 6500 cm~^. This double peak dif- 
fers from a linear combination of two-dimensional crystal 
and three-dimensional crystal peaks (see Fig. [3](c)) that 
would be expected naively, since both of its components 
appear at lower frequencies than the two-dimensional 
crystal and three-dimensional crystal peaks. While we 
have currently no interpretation for this intricate dou- 
ble peak, the observation of interface peaks calls for ad- 
vanced theoretical interpretations outside the scope of 
the present paper. In the long-wavelength limit below 
5000 cm^^, both the two-dimensional and inverse wood- 
pile photonic crystals are transparent. As a result, their 
reflectivities are low. Since the silicon flUing fraction is 
smaller for inverse woodpile photonic crystals than for 
two-dimensional crystals, they have a lower effective re- 
fractive index [43], resulting in a lower long- wavelength 
reflectivity for inverse woodpile photonic crystals. 

By combining this y-scan with six other scans for or- 
thogonal directions, crystal surfaces, and polarizations, 
we can measure the optical size of the crystal, in other 
words, the range over which the crystal has the same re- 
flectivity. This is in the x-dircction approximately 3 /im, 
in the y-direction 7 fj,m and in the z-direction 4 /im. The 
optical size is smaller in the x-direction than in the y- 
direction, which is in agreement with detailed structural 
information derived in Ref. |38j for inverse woodpile pho- 
tonic crystals. 



B. Overlapping stopbands for orthogonal 
polarizations and crystal surfaces. 

We have studied the position-dependent reflectivity 
spectra for two orthogonal crystal surfaces and two or- 
thogonal polarizations, by taking scans similar to the one 
presented in Fig. [3] From these scans we have selected 4 
(=2x2) spectra that are reproducible on several loca- 
tions on the crystal, see Fig. ^a-d), black curves. Fig. 
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Figure 3: (color) Position- and frequency-resolved reflectivity for an inverse woodpile photonic crystal, measured parallel with 
the rX-direction for _L-polarized light, (a) Microscope picture of the crystal placed in the setup. The bright spot near the 
center of the crystal is reflected focussed light of the supercontinuum white light source. The spot was moved over the crystal 
surface along the y-direction. (b) Reflectivity collected as a function of the y-position of the spot on the photonic crystal surface. 
The scan started on the surrounding two-dimensional crystal {y = /im), continued over the inverse woodpile photonic crystal 
(6 < y < 14 pm) and ended on the two-dimensional crystal. The spot was moved in steps of 0.46 ± 0.05 fim. The grey lines on 
the horizontal axis indicate the boundaries of the crystal. The white dashed lines indicate positions for representative spectra 
displayed in (c). (c) Reflectivity spectra measured at the y-coordinates given by the white dashed lines in (6). The calculated 
stop gap in the FX-direction for ^ = 0.19 is marked by the blue bar. 



|4|^a) shows the reflectivity centred in the FX-direction for 
_L-polarized light, and Fig. ^b) for || -polarized light, Fig. 

shows the reflectivity centred in the FZ-direction for 
_L-polarized light, and Fig. Qc?) for || -polarized light. A 
broad reflectivity peak appears for both directions and 
polarizations with maxima up to 67%. The stopbands 
overlap between 5900 and 6900 cm~^ (blue bar), corre- 
sponding to a relative bandwidth of 16%. This position- 
independent overlapping stopband for orthogonal polar- 
izations and crystal directions is a signature for the pho- 
tonic band gap with a relative bandwidth up to 16%. 

The reflectivity obtained from the FX-direction is 
shown in Fig. |4|a — 6). Below 5000 cm~^, the reflectivity 
behaves similarly for _L- and ||-polarized waves. Based 
on the filling fraction of this crystal, one would expect a 
reflectivity of approximately 10%, which is in very good 
agreement with our observations. A major reflectivity 
peak is observed for both polarizations; for _L-polarized 



light (Fig. ^a)) the peak appears between 5800 and 
7300 cm^^, for || -polarized waves the peak appears be- 
tween 5300 and 7100 cm~^. At higher frequencies the 
reflectivity varies around 20%. Reference spectra (grey) 
were collected on the silicon surface within 5 /zm distance 
from the inverse woodpile photonic crystal. The reflec- 
tivity varies between about 20 and 30%, which agrees 
closely with 31% reflectivity based on Fresnel reflection 
and a refractive index of n = 3.5 for a flat silicon sur- 
face. The reference spectra confirm that the observed 
stopbands in the inverse woodpile photonic crystal are 
indeed caused by the geometry of the nanostructure. 

The refiectivity obtained from the FZ-direction is 
shown in Fig. [4|c — d). Below 5000 cm~^, the reflectivity 
behaves similarly for _L- and ||-polarized waves. For the 
FZ-direction, fringes appear and the average reflectivity 
is higher than for the FAT-direction. These fringes are 
Fabry-Perot interferences caused by internal reflections 
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Figure 4: (color online) (a — d), Measured reflectivity (black 
lines) for an inverse woodpile photonic crystal along two or- 
thogonal directions (centred on FX and VZ) and orthogonal 
polarizations (perpendicular _L and parallel || to the pores 
perpendicular to the measurement axes). The overlapping 
frequencies of the stopbands (blue bar) reveal the presence of 
a photonic band gap. Reference reflectivity (grey lines) are 
measured on non-photonic parts of the sample within 5 fim 
from the photonic crystal, (e), Calculated bandstructure for 
both FX and TZ reveal a stopgap (light green bar), which is 
part of the photonic band gap (dark green bar). 
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Figure 5: (color online) Calculated stopgap for the VX- and 
rZ-directions as a function of relative pore radius ^. The 
band gap (light blue) is part of the broad stopgap (dark blue). 
The grey bar is the value of ^ based on parameters deter- 
mined from SEM images. The bandstructures in Fig. [2] and 
|4](e) were calculated at ^ = 0.19 (vertical dashed line) that 
agree best with the experiments. The horizontal dashed lines 
mark the edges of the common stopband determined from 
Fig. %a-d). 



in the inverse woodpile photonic crystal between the air- 
crystal boundary and the crystal-silicon boundary. The 
fringes are also present on the non-photonic reference 
spectra (grey) that are collected parallel to the first set of 
pores that form the surrounding two-dimensional crystal. 
Above 5000 cm~^ a major reflectivity peak is observed 
for both polarizations near the expected band gap. For 
_L-polarized light (Fig. |4]jc)) the peak appears between 
5900 and 6900 cm^^, for ||-polarized waves the peak ap- 
pears between 5300 and 7200 cm^^. The amplitude of 
the reflectivity peaks in the r.Z-direction are compara- 
ble with the reflectivity peaks of the FX-direction. For 
the FZ-direction, several troughs appear in the spectra, 
which are likely caused by surface contamination. 

If one overlays the spectra for _L-polarized light, com- 
pare Fig. |4]ja, c), it appears that the peaks overlap at the 
low-frequency side. In addition, for both the VX- and 
FZ-directions, the stopband is broader for ||-polarized 
waves than for _L-polarized waves. This is an experimen- 
tal indication that the pore radii of both sets of pores are 
identical, as expected from the symmetry of the struc- 
ture. 



V. DISCUSSION 

In order to interpret our observations, notably of the 
overlapping stopbands for two orthogonal crystal direc- 
tions and two orthogonal polarizations, we have per- 
formed photonic bandstructure calculations. In the 
course of the calculations, we have come to realize that 
pore radii determined from electron microscopy are an 
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overestimate of the true pore radii. This is reasonable, 
since parameters determined from electron microscopy 
have un uncertainty, typically near 2 to 3 % [IS]. There- 
fore, we have performed bandstructure calculations as a 
function of ^. Fig. |5] shows the calculated stopgap in the 
TX- and FZ-direction and the band gap versus ^. The 
broad stopgap corresponds to the intense stopbands in 
our reflectivity spectra. The horizontal dashed lines mark 
the edges of the overlapping stopband observed in Fig. 
l4l The grey bar is the value of ^ based on parameters 
aetermined from SEM images. For a relative pore radius 
between - = 0.19 and - = 0.20, the observed stopband is 
in agreement with the calculated band gap and stopgap, 
which is close to the value of - based on SEM images. 
This indicates that the crystal has an ^ between 0.19 and 
0.20. The corresponding bandstructure for ^ = 0.190 is 
shown in Fig. [5] and the bandstructure in the FX- and 
FZ-directions are shown in Fig. ^e). A stopgap is pre- 
dicted between 5727 and 7061 cm^^, which is part of the 
calculated band gap between 5727 and 6668 cm^^. The 
calculated band gap supports the observed signature of 
the photonic band gap between 5900 and 6900 cm~^. 

We have performed reflectivity measurements with an 
objective with NA= 0.65, resulting in a spectrum angle- 
averaged over OAin ± 10% sr solid angle in air. We have 
measured the reflectivity from two orthogonal crystal sur- 
faces and we observe intense peaks. Therefore, as argued 
above from the high reflectivity and the finite size of the 
photonic crystals, it is most likely that the observed stop- 
bands extend over 0.887r±10% sr solid angle. By invoking 
the symmetry of the crystal, the stopbands centred along 
the opposite —FX- and — FZ-directions are identical to 
the FX- and TZ directions. Therefore, our observations 
would correspond to stopbands up to 1.767r± 10% sr ex- 
ternal solid angle of the inverse woodpile photonic crys- 
tal, in other words nearly half of all directions. In previ- 
ous experiments on inverse opals without bandgap !46] , 
nearly half of all directions were found to be excluded 
for propagation for a vanishingly narrow frequency band. 
Therefore, the present result is to the best of our knowl- 
edge the largest solid angle for which a broad photonic 
stopband has been reported. 

The observed common stopband between 5900 and 
6900 cm~^ is well explained with bandstructure calcu- 
lations. A limitation to our interpretation is that we 
invoke a theory for infinite crystals. Theoretical cal- 
culations of reflectivity for finite crystals could result 
into new insights in the obtained spectra, for instance 
how light propagates inside photonic band gap crystals, 
which is an extremely complex problem as described in 
Ref. [37]. A well known method for finite crystals are 
finite-difference time domain simulations |48| . Unfortu- 
nately, however, full three-dimensional simulations with 



reasonable parameters for our experiments will require 
very extensive numerical calculations [49] that are out- 
side the scope of our present experimental study. Scat- 
tering matrix methods [501 ^FT provide a good descrip- 
tion of three-dimensional crystals with finite thickness. 
However, the method assumes that a crystal is infinitely 
extended in both transverse dimensions which is not the 
case for our crystals. Hence we propose that a theory for 
three-dimensional finite crystals is called for. It will be 
even more difficult to implement the experimental bound- 
ary conditions into such calculations, such as the shape 
of the pores and the NA of the microscope objective. At 
any rate, the agreement between our calculations and our 
experimental data is gratifying. 

For good measure, the observed signature of the pho- 
tonic band gap is not a complete experimental observa- 
tion of such a gap. At this time, we can not experi- 
mentally exclude that there are crystal directions where 
the stopbands are either closed or shifted to different fre- 
quencies. While we consider these possibilities unlikely 
based on the theoretical study of our crystals and in view 
of the measured crystal parameters, a fully experimental 
observation of a photonic band gap must entail an exper- 
imental method where the local density of optical states 
is probed. An example would be the doping of the crys- 
tals with quantum dots 52^ and measuring the emission 
rates [SHE]- 



VI. CONCLUSIONS 

In conclusion, we have studied the polarization- 
dependent refiectivity centred on the FX- and TZ- 
direction of three-dimensional silicon inverse wood- 
pile photonic crystals. We have observed a position- 
independent overlapping stopband for different polariza- 
tions and directions, which are supported with calculated 
bandstructures. This is the signature of a photonic band 
gap with a bandwidth up to 16%. 
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